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Abstract 

In the present paper, we carry out a detailed analysis of the presence and mixing of various 
families of collective bands in 188 Pb. Making use of the interacting boson model, we construct 
a particular intermediate basis that can be associated with the unperturbed bands used in more 
phenomenological studies. We use the E2 decay to construct a set of collective bands and discuss 
in detail the B(E2)-values. We also perform an analysis of these theoretical results (Q, B(E2)) 
to deduce an intrinsic quadrupole moment and the associated quadrupole deformation parameter, 
using an axially deformed rotor model. 

PACS numbers: 21.60.Ev:21.60Fw 
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I. INTRODUCTION 



The lead isotopesprovide a unique laboratory to study the phenomenon of shape coexistence 



3 ,|4| and have v ery recently been subject of much experimental and theoretical 

The combined effect of the proton shell- 



in nuclei ll|, z , 
interest p, |6 , 

gap at Z=82 for the Pb-nuclei and the large number of valence nucleons outside the closed 
iV=126 core (in this case, neutron holes), results in an important lowering of the energy of 
proton particle-hole excitations Q|. More specifically, near neutron midshell (at A=104), 
the proton 2p-2h and 4p-4h excitations descend to very low excitation energy because of 
the very large proton-neutron binding energy that results from the interactions between 
the core-excited protons across the Z=82 shell closure and the large number of valence 
neutrons outside of the A=126 shell closure. As a consequence, mixing can result between 
various families of configurations having approximately the same excitation energy More in 
particular, the lowest-lying + and 2 + states can become strongly mixed such that it is very 
difficult to assign a 'configuration label' to them. In section |HJ we succintly describe the 



essentials of the configuration mixing using the interacting boson model (IBM) [I8l.ll9l.l20j] .In 
section 11111 we have in mind to obtain a better understanding in the nature of the low- 
spin collective states and to construct particular bands by making use of the calculated E2 
decay probabilities starting at the high-spin states in the particular case of 188 Pb. In a final 
section|lVj we use the calculated values (using the IBM) for quadrupole moments and B(E2)- 
values in these bands in order to extract collective model parameters (intrinsic quadrupole 
moments, quadrupole deformation ft Q ) in order to highlight the equivalence between the IBM 



approach, used as a highly-truncated shell-model calcu 
model, as can be derived from mean-field methods ((J |s , 



ation, and a geometrical rotational 







16, 



21, 



221. 



II. CONFIGURATION MIXING 



In a recent study that concentrated on describing intruder bands and configuration mixing 
in the neutron deficient Pb-isotopes, a three-configuration mixing calculation has been per- 

23j . We refer to that paper for 



formed in the context of the interacting boson model (IBM) 
more details but present succintly the main points. Approximating the regular and intruder 
states as built from proton Op-Oh, 2p-2h and 4p-4h excitations across the Z=82 proton closed 
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shell mainly that interact with the large number of valence neutrons outside of the N=126 
neutron closed shell, the Hamiltonian within the IBM takes the form 



with 



H — H reg + H 2p -2h + #4p-4/j + Vmix , (1) 
H r eg ~regTld ~\~ l^regQreg ' Qreg ; (2) 

Hi = Eih d + KiQi -Qi + Ai, (3) 



and 



Vmix Vmix,l Vmix,2 j (4) 

V mi x,i = oa (s f • s f + s ■ s) + (3i (d) ■ d) + d ■ d) . (5) 
The quadrupole operator Qi is defined as 

/ ~ \ (2) / ~\ (2) 

Qi = (sU + d^s) + Xi (d ] d\ . (6) 



For details on the notation we refer to Fossion et al. |23J. The diagonalisation of the energy 
matrices, corresponding to the Hamiltonian (fljl. is carried out in the t/(5)-basis, expressing 
the eigenvectors in the [iV]©[iV" + 2]©[iV" + 4] model space. Unfortunately, using this method, 
one has no clear insight in the amount of mixing between the unperturbed bands that 
result from a diagonalisation in the separate subspaces [TV], [iV+2], and [iV + 4], respectively. 



So instead of a complete diagonalisation of the Hamiltonian matrix, we first rotate 
into an intermediate basis in which the separate parts of the Hamiltonian i.e. H reg , 
H2 P -2h, and H^-^h only, become diagonal. We thus have three different bases: 



• The U(5)-basis |J, k)^, with J the angular momentum, N the number of bosons and 
k the rank number, 

• the basis in which the full Hamiltonian (0) is diagonal | J, i) , and, 

• the intermediate bases in which the Hamiltonian (excluding V m i X ) is diagonal in 
the three different subspaces [N] , [N + 2] , and [N + 4] . Respectively they are denoted 
as \J,l)' N , \ J,l)' N+2 , and \J,l)' N+4 , with I also a rank number. 
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The above bases are connected in the following way. The basis in which the full Hamiltonian 
is diagonal (expressed in the C/(5)-basis) reads 

\J,i)= E a £( J )l J '^ 

k=l 

dim^f+dimffj_2 

+ E <t\j)\j,i) N+2 

/=dimjv+l 
dim jv +dimjv + 2 +<i* m JV+4 

+ E a^ 4 (J)|J,m) w+4 , (7) 

m=dimjv+rfiT»jv+2+l 

and the 'intermediate' basis becomes 

\j,iy N = EWi J .^- ( 8 ) 
fc=i 

(similarly for iV + 2 and iV + 4). Dim^, dim^r+2, and dini7v+4 are the dimensions of the 
corresponding configuration spaces for a certain angular momentum J containing N, N + 2, 
and iV + 4 bosons respectively. The matrix B N diagonalises the configuration with N bosons 
and the matrix A diagonalises the full Hamiltonian (0) directly (the indices N, N + 2 ,and 
iV + 4 for matrix A are only added for the sake of clearness). We omit the dimensions of 
the summations from now on. 

Rotation of the Hamiltonian matrix expressed in the ?7(5)-basis into the 'intermediate' 
basis results in the energy levels of a set of bands in the Op-Oh, 2p-2h and 4p-4h subspaces 
separately (see also left part in Fig. Q). These bands correspond to the unperturbed bands 
that are extracted in phenomenological calculations as carried out by Dracoulis et a/.j^|, 
Allatt et al. j^j], and Page et al. j^j]. From now on we will call the energy levels (bands) 
resulting from rotation of the Hamiltonian matrix into the 'intermediate' basis unperturbed 
levels (bands) in order to avoid confusion and we will denote them as the | J, L)' N states (see 
eq. ®). 

By calculating the mixing matrix elements of V m i X in this 'intermediate' basis, we obtain 
the mixing matrix elements for all unperturbed levels. The knowledge of these unperturbed 
bands and their mixing matrix elements makes the process of configuration mixing in a 
nucleus more transparent than before. Starting from the experimental level energies at high 
spin and the corresponding moment of inertia, a set of unperturbed experimental bands 
can be deduced extrapolating to the low-spin members of these bands 



xperimental banas 



«i = fa a 2 = fa Ai 



8.5 keV 23.4 keV 1923 keV 



TABLE I: IBM-parameters used for 188 Pb. 



The comparison between the IBM unperturbed bands and these extrapolated unperturbed 
experimental bands, together with the knowledge of the full energy spectrum and the 
B(E2)-values form an extensive test for the parameters that describe a certain isotope chain. 

Starting from IBM-parameters for the Pb-isotopes as determined by Fossion et al. 
(23!, a slightly different fit was performed. The parameters for H reg , H 2p -2h, and _ff 4p _ 4 /j 
remain unchanged, except for e reg which was taken 0.92 MeV instead of 0.90 MeV. The 
mixing parameters and fa were fixed in 196 Pb. The value for Ai was obtained as the 
result of a fitting procedure for the Pb-isotopes (A— 186- 196). Basically Ai was fitted for 
186 Pb and 196 Pb and the Ai for the other isotopes was chosen following a linear variation 
between Ai( 186 Pb) and Ai( 196 Pb). Then, the difference between the experimental 2 and 
the IBM 0^~ was taken for all isotopes considered and added to the corresponding Ai. This 
method gives a better description of the slope of the energy levels through the isotope chain. 
Moreover, the Ai( 188 Pb) obtained in this way is in good agreement with the theoretical 
prediction that makes use of experimental separation energies Q]. Since A 2 is associated 
with the unperturbed energy to excite 4p-4h configurations, the value was taken as 2-Ai 
j3| . The parameters for 188 Pb are listed in Table |1J The low-energy part of the resulting 
IBM-spectrum for 188 Pb is presented in Figure [0 The comparison with the experiment will 
be discussed in the next section. 



By inspecting the unperturbed lowest bands (see left part of Fig. |TJ) one notices 
that mixing modifies the nature of the 2p-2h and 4p-4h + states and the three 2 + states 
because of the small energy differences between these unperturbed states (see Figure [IJ. 
The interaction matrix elements between the first three + states and the first three 2 + 
states are given in Table |Hj We now make the approximation to neglect the mixing between 
the lowest and first excited + states. We can make this approximation because the mixing 
matrix element coupling those two states is less than half the magnitude of the interaction 
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FIG. 1: Low-energy part of the IBM-spectrum for Pb. The three configurations taken from the 
left show the absolute energies for the lowest unperturbed bands. At the extreme right, the IBM 
spectrum after diagonalisation of the full Hamiltonian Q is shown. 
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TABLE II: Mixing matrix elements for the first three unperturbed + and 2 + states. The matrix 
elements are expressed in MeV. 



matrix element coupling the Of and 0^" states. This is strengthened by the fact that the 
energy difference between the Of and 0^ states is much larger than between the and 
Ot states. The situation thus reduces to a simple two level mixing problem for which the 

n 

energies and eigenvectors are easily determined j29(. We obtain a level at 494 keV and one 
at 882 keV, compared to the absolute energies 489 keV and 853 keV resulting from the 
full diagonalisation. Although we neglected the effect of the other 55 + levels, we obtain 
a good result for the energies using two-level mixing only. This is also reflected in the 
wavef unctions. According to the two-level mixing problem, the wavefunctions look like 

ipi = an0i O) + a2i02° } > ( 9 ) 
^2 = a 12 <ri 0) + a 22 # , (10) 
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with 

an = =R22 ,012 = ±021 • (11) 

Table IIHI gives the coefficients for the wavefunctions of the + (and 2 + ) states resulting 
from the full diagonalisation expressed in the unperturbed basis (jHJ)- Only the coefficients 
of the lowest three unperturbed + (and 2 + ) states are shown. Inspecting the coefficients 
of the unperturbed 2p-2h and 4p-4h + wavefunctions \0f)' N+2 and 10^)^+4 respectively 
for the 0^ and O3" wavefunctions, we find the same structure for the magnitudes as in 
equation One notices that a simple two- level mixing approximation cannot reproduce 
the phases resulting from the full diagonalisation. Calculation of the coefficients in the 
two-level approximation for the + states results in the values 0.7696 and 0.6384, which are 
in good agreement with the magnitudes resulting from the full diagonalisation. 





10^ 
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1^1 }'n+4 


19+V 


\^iYn+2 


|2i }'n+4 


|of> 


-0.9897 


-0.1310 


0.0309 


|2f> 


0.1907 


-0.6379 


0.7091 


\ot) 


0.1216 


-0.7624 


0.6148 




0.5964 


-0.4906 


-0.5664 


|o 3 + ) 


-0.0611 


0.5893 


0.7655 




0.7696 


0.5114 


0.2187 



TABLE III: Coefficients for the wavefunctions of the + (and 2 + ) states resulting from the full 
diagonalisation expressed in the unperturbed basis (jHJ). Only the coefficients of the lowest three 
unperturbed states \0f)' u and (2^)^ (v = N, N + 2, N + 4) are given. 

The above results indicate that the two-level approximation to the full diagonalisa- 
tion is a very good one. It thus follows that a knowledge of the mixing matrix elements and 
the energies of the unperturbed levels can lead to a better understanding in the process of 
configuration mixing. We like to point out that this severe approximation cannot reproduce 
the correct phase, so we cannot use these approximated wavefunctions for the calculations 
of B(E2)-values. 

III. B(E2)-VALUES AND CONSTRUCTION OF COLLECTIVE BANDS 

In the process of separate diagonalisation of the full Hamiltonian (0) (excluding V m i X ) in 
the three subspaces separately, we have also made an 'intermediate' calculation concerning 
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B(E2)-values. The knowledge of the basis states (J7J) and (JHJ) gives rise to two interesting 
expressions for the reduced matrix element of a transition between the initial Ji(i) and final 
Jf{f) state, i.e. 



(J f , f 11^2)11^,1) = 



k p 



=N 



; q 



u=N+2 



+ E E <i( J *KA J f) AJf,r\\T(E2) \\J h m) v 



u=N+4 



(12) 



and, 



(J,,f\\T(E2)U,i) = E [,<J/y||T(S2)||J 4 , a >i 



k,p,s,s' 



l,q,t,t' 



u=N 



v=N+2 



m,r,u,W 



u=N+4 

(13) 

with component of transposed matrix B of equation (jHJ). Expression (JT3j) is most 
interesting because it allows us to check which transitions in the unperturbed bands make 
up for an important contribution to a certain transition Ji(i) — > <//(/). We now apply the 
above method to the particular case of 188 Pb. 



For 188 Pb, only two experimental B(E2)-values are known 11 1 

B(E2;2t -> 0+) = 5(3) W.u. , 
E(£2;4+ -> 2+) = 160(80) W.u. 



In the calculation of the E2-transition rates, we use the consistent-Q procedure |30|] to 
determine the E2-transition operator as 
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(14) 
(15) 



T(E2) = E ej (s f d + d f s) (2) + Xi (t^d 
i=i L 



~\(2) 



(16) 



Xreg X2p-2h X4p-4h &reg &2p-2h — ^4p-Ah 

0.515 -0.680 0.11 0.132 

TABLE IV: Parameters for E2-transitions in 188 Pb. The effective charges are expressed in e • b, the 
X are dimensionless. 

So we choose the values for X2 P -2h and Xip-ih as obtained in j^J and fit the effective charge 
to those two known data. We took t<ip-2h and e^p-Ah 1-2 times e reg . Table ITVl lists the 
parameters used. We point out that we have chosen Xreg = 0. We made this choice because 
there are no data available for transitions within the regular band, hence we are not able 
to fit Xreg to known experimental values. Evidently, the groundstate is regular and we 
know the B(E2)- value for the 2+ — > Of transition, but as Xreg is of no influence to 2 + — > + 
transitions in the [7(5)-limit, we cannot make use of this transition in determining Xreg- 

For this choice of parameters the IBM calculation yields the following results 

B(E2; 2+ -> g.s.) = 0.0195 e 2 b 2 or 3 W.u. , (17) 
B(E2; 4+ -> 2+) = 0.9747 e 2 b 2 or 152 W.u. . (18) 



The fact that the parameters are able to reproduce two B(E2)-values that differ two orders 
of magnitude, permits us to make further theoretical predictions. 



To begin with, it is interesting to make a comparison between the B(E2)-values for 
the unperturbed 2p-2h and 4p-4h band and the theoretical estimate using a pure rotational 



band 



31 



32] 



B(E2] Jf + 2 g.s.band -> J f g.s.band) 



/ 

J f + 2 2 J f 


3 (J / + l)(J / + 2) 



c (2i > + 1; 



c 



(19) 



(20) 



2(2J / + 3)(2J / + 5) ' 

with C a constant factor depending on the nucleus. Since C is unknown, we can only 
compare with the B(E2)- values relative to a certain transition. In Figure 121 we illustrate 
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FIG. 2: The ratio B(E2;J f + 2 -> J f )/B(E2; 10+ 8+) in the unperturbed 2p-2h and 4p-4h 
bands compared to the theoretical estimate for a pure rotational band. The full line denotes the 
theoretical estimate, the dot dashed line the ratio for the 2p-2h unperturbed band and the dotted 
line the ratio in the 4p-4h unperturbed band. 

the .//-dependent part of equation (}2T)j) and draw a parallel with the B(E2)-values in the 
unperturbed bands, all relative to the 10 + — > 8 + E2-transition. 

One can clearly see that the band built upon the unperturbed 4p-4h + state follows 
the rotational structure very well, while the band constructed upon the unperturbed 2p-2h 
+ exhibits larger deviations, especially for the transitions between the lower spin states. 
Because x^ P -ih was fixed using I-spin symmetry j^J, we can explain the rotational- like 
behaviour of the unperturbed 4p-4h band by comparing with the re gula r band in 180 W. 
This band is recognised as a K n = + rotational groundstate band 33(, which explains 
the good agreement of the slope of the ratio B(E2; Jf + 2 — > Jf)/B(E2\ 10 + — ► 8 + ) in the 
unperturbed 4p-4h band with the theoretical one. 



Finally, we come to the question about labelling the mixed states, resulting from 
diagonalisation of the full Hamiltonian (|T|). in a meaningful way into a given band. As the 
criterium, we start from the calculated high-spin members of the two bands and follow the 
E2-decay. We define the collective band through a sequence of E2-transitions down to low 
spin such that the intra-band B(E2)-values are bigger than the inter-band B(E2)-values. 
Moreover in the case of 188 Pb, we have as a second criterium the fact that the intra-band 
transitions between higher spin (and thus relatively unperturbed) states should exhibit 
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transition 


B(E2)-value 


4+ -» 2+ 
^2 ^ z 3 


44 W.u. 


2+ -» 2+ 


36 W.u. 


^3" — * ^2 


2 W.u. 



TABLE V: Inter-band transitions involving the 2g state 




FIG. 3: Experimental and theoretical level scheme in 188 Pb. The arrows denote the B(E2)-values 
for a given transition, expressed in W.u. The experimental data were taken from 
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13 



34] 



an approximately constant behaviour as illustrated in Fig. |21 The final results are shown 
in Figure EH In Table all inter-band transitions involving the IBM 2^ at 969 keV (not 
drawn in Fig. |3J) are shown. The transitions between Band II and Band III in Figure |3] 
that are not shown, are less than 1 W.u. We like to stress that the 725 keV + state in 
Figure H3 is indicated with 



two experimental studies 



24, 



shed line and tentative energy. This state was identified in 
but recently performed experiments QQ do not reach 
conclusive evidence for this state and its energy. The position and existence of this level 
thus is still under discussion. 
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Comparing the IBM results with the experimental bands, we notice that the struc- 
ture of the collective bands is reproduced rather well. The energy levels that constitute the 
more collective Band II are in good agreement with the experimental results. The structure 
of the Band III is also described rather well although the value of A2 P -2h seems a bit too 
small and the mixing parameter between the two intruder configurations slightly too large. 
Although the 0^ and the Ojj" state are placed in the collective bands the other way round in 
reference [ll]], here we switched them. This change was performed following the definition 
that inter-band transitions must be small compared to intra-band transitions. In the other 
case than the one depicted here, the IBM-calculations give transitions from the collective 
Band III to the regular ground state that are larger and thus more collective than the 
transition to the + state of Band III itself. Such a decay pattern should definitely be 
discarded. The only way to obtain a consistent picture in which the E2-transitions exhibit 
a collective nature within the collective bands and give rise to small E2-transitions between 
the different bands, is by associating the 0^ level to Band II and the Ojj" level to Band III. 

When we take a closer look at the E2-transitions between the low-lying levels, some- 
thing interesting occurs. Using expression (|13|). one can single out those contributions of 
the reduced matrix elements for the transitions in the unperturbed bands (thus reduced 
matrix element and coefficient!), that make up the major contribution to a given transition 
Ji(i) — > J/{f) with the condition that their sum may not deviate more than 10% from 
the total value with all contributions taken into account. In the case of 188 Pb, it turns 
out that for the so constructed Band II and Band III, the main contribution to intra-band 
transitions always consists of the term with the corresponding lowest rank transition in the 
unperturbed 2p-2h band and the one in the unperturbed 4p-4h band. Thus 

(J f ,f\\T(E2)\\J i7 i)<* 

l,q 

m,r 

= Z(N + 2)(J f ,f\\T{E2)\\J i ,i) + Z(N + 4:)(J f J\\T(E2)\\J i ,i) (21) 
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FIG. 4: (colour online) The left part shows the ratio's Z(N + 2) and Z(N + 4) (see equation l|21|)) 
for the intra-band transitions in Band III, the right part shows the ratio Z(N + 2) and Z(N + 4) 
for the intra-band transition in Band II. 

with \Ji,i) and \Jf,f) both states in Band II or both states in Band III. This fact is as 
expected. More interesting is the ratio Z of those two main contributions to the total value 
of the reduced matrix element. These results are depicted in Figure |U 

Notice that the ratio Z(N + 4) becomes larger than the ratio Z(N + 2) for the 
2\ — > 0^ transition in Band III. This seems to be in contradiction with the fact that we 
placed the and Ojj" states in Band III as for all other transitions in this bands Z(N + 2) 
dominates over Z(N + 4). However, all other orderings of these levels into the bands (II 
and III) are in contradiction with the E2-decay pattern and our definition of intra-band and 
inter-band transitions strengths in order to define a collective band. This very large ratio 
Z(N + 4) for the 2^~ — > 0^ transition in Band III is due to the combined effect of the large 
mixing and the relative magnitudes of the two contributing reduced matrix elements (j2"Tj) . 

So we can conclude that the only physically meaningful way to characterise the low- 
lying states is by constructing the different bands using the intensities of the sequences of 
E2-transitions. Attempts to label the extremely mixed low-lying states as '2p-2h' or '4p-4h' 
lead to inconsistencies. 
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IV. COMPARISON WITH THE COLLECTIVE ROTATIONAL MODEL 



Having constructed two different collective bands starting from the IBM and using the 
prescription to define bands on the basis of the calculated B(E2)-values, it is very useful 
to check for consistency with the results of other theoretical approaches. First of all we 
concentrate on the quadrupole moments. 

Within the IBM, the quadrupole moments are calculated as 



Q(J) = \I^ L (JJ\T(E2) \JJ) 



16tt / J 2 J 



(J\\T(E2)\\J) . (22) 



5 I -J J , 

From the point of view of the collective rotational modelj^. I32I Is^ . the electric quadrupole 
moment is defined as 



MRtt 

Q = \/ —(J,K,M = J\M(E2,0)\J,K,M = J) . (23) 



For K = bands this reduces to 



W " (J+l)(2J + 3) Wo ' [ZV 



with the intrinsic quadrupole moment Qq defined as 



Ql = \J l -^{K = Q\M{E2MK = 0) . (25) 

Equating the quadrupole moments of equation (|23*|) and allows us to extract an equiv- 
alent intrinsic quadrupole moment Qq. Likewise one can use the IBM B(E2)-values to 
extract an equivalent intrinsic quadrupole moment Qq using the collective rotational model 
B(E2)-expression 

B(E2; ( J + 2), K = -> J, K = 0) = 

5 / J + 2 2 J \ 

(2J+1) (Q°oY- (26) 



167T 







Having extracted an intrinsic quadrupole moment, one can deduce a deformation parameter 
(3q using the expression 

Qo = 3 4=P° ■ ( 2? ) 

V 57T 
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If the concept of a collective band characterised by a single Qq value is to make sense, 
one expects that values of Q® extracted with these two procedures will not differ much. 
Moreover, one expects only a moderate variation of Q$ as a function of J along the band. 




| -6 ■ 

a 

-8 I 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 

J 

FIG. 5: Intrinsic quadrupole moment extracted from transition matrix elements, calculated within 
the IBM. The intrinsic quadrupole moment is expressed in e • b. Results for the 2p-2h unperturbed 
band are represented with •, results for the 4p-4h unperturbed band with ■. Bands II and III are 
represented respectively by * and ▲. 

In Figure El we have plotted the intrinsic quadrupole moment Q® derived from the 
E2 diagonal matrix elements. We clearly see that the intrinsic quadrupole moment stays 
approximately constant for the unperturbed bands. The intrinsic quadrupole moments of 
the unperturbed 4p-4h band are positive, indicating a prolate deformation of the nucleus, 
while the negative sign of the intrinsic quadrupole moments in the unperturbed 2p-2h band 
is consistent with an oblate deformation. The values of the intrinsic quadrupole moments 
for the unperturbed bands obtained from the IBM diagonal transition matrix elements 
are in good agreement with the quadrupole moment magnitudes reported by Dracoulis 
et al. [9|. They calculated a magnitude of 3.4 eb for the intrinsic quadrupole moment 
of the unperturbed oblate band, while the IBM gives an intrinsic quadrupole moment 
varying between -3.4 eb and -3.7 eb. The same is valid for the unperturbed prolate band 
with a quadrupole moment of 6.2 eb resulting from their band mixing and a quadrupole 
moment varying between 5.6 eb and 6.6 eb resulting from the IBM-calculations. The 
mixing causes the bands to become less prolate and oblate so the intrinsic quadrupole 
moments for Band II and III respectively are considerably smaller, as can be seen in Figure^ 
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Extracting intrinsic quadrupole moments from B(E2)- values gives a similar picture 




2 4 6 8 10 12 14 16 2 4 6 8 10 12 14 16 



J J 

FIG. 6: The upper figure shows the comparison between the magnitudes of the intrinsic quadrupole 
moments extracted from the theoretical B(E2)-values (•) and the magnitudes of the intrinsic 
quadrupole moments extracted from the diagonal matrix elements (■), both for the unperturbed 
2p-2h band. The lower figure shows the same comparison for the unperturbed 4p-4h band. 

as in Fig. El More interesting is the comparison between the intrinsic quadrupole moments 
calculated from the B(E2)-values and the absolute values of the intrinsic quadrupole 
moments as extracted from equations (|23|) and (|24|) . Figure El depicts the comparison for 
the unperturbed 2p-2h (upper part) and the unperturbed 4p-4h band (lower part). It can 
be clearly seen that the unperturbed 4p-4h band exhibits a rotational-like behaviour to 
large extent, while the unperturbed 2p-2h band is less rotational. This is in agreement with 
our conclusion for Fig. |21 

Finally, we have also extracted the deformation parameter (3q. For the unperturbed 2p-2h 
band /3q varies between -0.11 and -0.12, while for the unperturbed 4p-4h band (3q stays 
between 0.19 and 0.22. These values are smaller than the deformation parameters calculated 
by Bender et al. p. 

V. CONCLUSION 

In the present paper we have carried out a detailed analysis of the presence of various 
families of collective bands, in particular for 188 Pb. We have started from an algebraic model 
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approach (the interacting boson model) in order to truncate the extended shell-model space 
that incorporates the presence of proton particle-hole excitations across the Z=82 closed 
shell. Moreover, making use of a concept called intruder symmetry, we have been able to 
reduce the number of parameters in the present description. A first calculation has been 
carried out by Fossion et al. [23| , accentuating the presence of three different families. Here, 
we have defined an intermediate basis that defines three separate systems by diagonalizing 
the Hamiltonian in the Op-Oh, 2p-2h, and 4p-4h configuration spaces separately. This basis 
allows to understand the mixing between the unperturbed bands in a transparent way. 
Moreover, we have used the E2-decay, starting at high spin, to define "physical" bands also 
progressing to low-spin members. Here the conclusion points towards an important mixing 
between the 0^~, O3" and the 2f, 2\ and 2% band members, still allowing the separation 
into two collective band structures. A simple reanalysis of these two bands within the 
collective rotational model is consistent with prolate and oblate band characteristics for 
the unperturbed 4p-4h and 2p-2h bands respectively. Extracted magnitudes of the intrinsic 
quadrupole and the deformation parameters are in good agreement with calculations 
starting from mean-field approaches. 

It is our aim to carry out a similar analysis for the other Pb-nuclei near the neutron 
mid-shell region and follow the mixing patterns between the two bands when moving 
away from the mid-shell region (at N=104) and to incorporate an extensive study of 
electromagnetic properties (E2- and EO-transitions) and the study of isotopic and isomeric 
shifts as well. 
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